The periodic Anderson model with correlated conduction electrons 
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We investigate a periodic Anderson model with interacting conduction electrons which are de- 
scribed by a Hubbard-type interaction of strength U c . Within dynamical mean- field theory the 
total Hamiltonian is mapped onto an impurity model, which is solved by an extended non-crossing 
approximation. We consider the particle-hole symmetric case at half-filling. Similar to the case 
U c = 0, the low-energy behavior of the conduction electrons at high temperatures is essentially 
unaffected by the / electrons and for small U c a quasiparticle peak corresponding to the Hubbard 
model evolves first. These quasiparticles screen the / moments when the temperature is reduced 
further, and the system turns into an insulator with a tiny gap and flat bands. The formation of 
the quasiparticle peak is impeded by increasing either U c or the c-f hybridization. Nevertheless 
almost dispersionless bands emerge at low temperature with an increased gap, even in the case of 
initially insulating host electrons. The size of the gap in the one-particle spectral density at low 
temperatures provides an estimate for the low-energy scale and increases as U c increases. 
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I. INTRODUCTION 

The usual explanation for the formation of heavy 
fermions in compounds with rare-earth or actinide ele- 
ments is based on the Kondo effect B u Thereby, the char- 
acteristic low-energy scale arises from the spin-screening 
of the local moments by a non-interacting electron gas. 
The periodic Anderson model is considered as the most 
promising candidate to at least qualitatively describe 
the rich physics of these materials. This standard sce- 
nario, however, fails to explain the heavy-fermion behav- 
ior found in the electron-doped cuprate Nd2- x Ce x Cu04 
discovered a few years ago.H In particular, the estk 
mated low-energy scale is orders of magnitude too small.cl 
Since undoped NdoCu04 is an antiferromagnetic charge- 
transfer insulatorQ'BL despite of one hole per unit cell, it 
has been suggested^ that this discrepancy is due to the 
strong interactions among the electrons introduced by 
doping. 

As regards the influence of correlated conduction elec- 
trons on the Kondo effect, up to now attention has been 
focused on the case of a magnetic impurity embedded in 
a correlated host which is either described by a Luttinger 
liquid in one dimensionLTu oiJpv-some kind of Hubbard 
model in higher dimensions OTJ In all of these cases 
one finds a strong dependence of the low-energy scale on 
the interaction strength of the conduction electrons and 
its increase with increasing interaction strength. For a 
lattice of moments hybridizing with correlated electrons 
only few results exist. taO 

As a first step towards understanding the effect of con- 
duction electron interactions on the formation of heavy 
fermions, we consider a lattice of /-electrons that hy- 
bridize with conduction electrons which themselves are 
correlated. These correlations will be described by a 



Hubbard-type interaction. The resulting model combines 
a periodic Anderson model with a Hubbard model. A 
limit in which one may obtain sensible results for this lo- 
cally highlyuicpjelated model is the limit of large spatial 
dimensions. EjpTJ In this limit the dynamics becomes es- 
sentially local.Ej Hence for any correlated model, a single 
(correlated) site may be chosen and embedded in an effec- 
tive medium which has to be determined self-consistently 
("dynamical mean-field theory!!) j— The model reduces to 
an Anderson impurity model.tJtj InfJaetp-iBJias been 
shown that besides the Hubbard mode£IE.al23'E2 the peri- 
odic Anderson model with uncocrelated conduction elec- 
trons is amenable to this limit .E3c2l 

In the next section we introduce the model and de- 
rive the corresponding impurity model. The impurity 
model is solvcd_niimerically by an extended non-crossing 
approximatiorE3~E3 which is derived in Sec. [njas well. In 
Sec. Ill we present results for the particle-hole symmet- 
ric case at half filling. Assuming a paramagnetic ground 
state we study the influence of weak correlations on the 
one-particle spectral density and discuss how the heavy 
bands emerge. At low temperatures a gap forms which 
in the free case is related to the Kondo temperatures 
and we discuss how the size of the gap depends on the 
strength of the correlations of the conduction electrons. 
We finally conclude in Sec. IV. 



II. MODEL AND METHOD 

A. The lattice Hamiltonian 

In the following we consider the simplest version of 
the periodic Anderson model and allow for interacting 
conduction electrons, 
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H = H C + H f + H cf (1) 
H c = ^(e fc + e c - M)4, CT c fc, ff + U c E < T <1 ( 2 ) 

fee - i 

»/ = (e/ - m) E /i,«r/i,a • '>E n Wi 

Here, /£P destroys (creates) an electron in the localized / 
orbital at site i with spin a, and Uf is the Hubbard inter- 
action of the localized / states (e/ < 0). The c operators 
refer to the conduction electrons which are described by 
a Hubbard model with an interaction U c being typically 
smaller than Uf. /i denotes the chemical potential and 
V measures the mixing between the c and / subsystems. 
In the following we will refer to the Hamiltonian ([!]) as 
"(periodic) Anderson-Hubbard model." 

In what follows we will concentrate on the one-particle 
spectra in a paramagnetic phase. Information on the 
underlying lattice will enter the dynamical mean-field 
equations onlv-via the density of states of the conduc- 
tion electrons .EJ For simplicity, we therefore consider a 
semicircular density of conduction electrons states with 
width 2D in Eq. (0) 



P(z) 



nD 



Z \2 



(5) 



which arises from hopping on a Bethe-lattice of coordi- 
nation number Z with matrix element t = D/(2y/Z) in 
the limit Z — > oo. We will use D = 1 as unit of energy 
throughout this paper. 

Given the non-interacting Green's functions, Go(k,z), 



Go X (M) 



-( z ~ ( e / ~ M) 
-V 



-V 

z - (ejfe - n) 



(0) 



and the full Green's function G(k, z), the self-energies are 
defined by Dyson's equation 



G~ 1 {k,z) = G ( 7 1 (fc,z)-E(fc,z) 



E(k,z) = 



_ ( E f (k,z) Ef c (k,z) 



E c f(k,z) S c (fc,z) 



(7) 
(8) 



B. The impurity model 

The dynamical mean-field theory assumes that the self- 
energy is a local quaoiita: wbirVi is correct in the limit of 
infinite dimensions .EJ'Ej'EScj The lattice model can then 
be mapped onto an impurity modelE3 which is seen as 
follows: The self energy (||) of the lattice model (|l|) is 
given as the derivative of a functional 3>[G] of the full 
Green's functions a 



^jj,v\iij\ z ) — 



(9) 



Here, i and j are real space coordinates and /z, v corre- 
spond to c and /. If the self-energy of the lattice model is 
local, $ depends on the local Green's functions, G(i,i), 
only. Thus, <E> can be generated from an impurity model. 
Solving the impurity model we know $ and, hence, the 
self-energy as functional of the impurity Green's func- 
tion Q: S[0] = 5Q/5Q. We now identify Q with the local 
Green's function of the lattice, 



G\ oc {z) 



jrT, G ^ z ^ I dep(e)G(e,z) . 



(10) 



Note that the fc-dependence enters only via e& into 
G(k,z), thus G\ oc (z) can be expressed by an energy in- 
tegration. From Gioc = G, we find the actual value of £ 
from the functional equation 



G loc (z) = / dep(e)G(e,z) 

dep(e)(Go\e,z)- S[G loc ]) _1 . (11) 



Technically, Eq. (11) determines the free Green's function 



of the impurity model which is of course not fixed by 

The Hamiltonian of the impurity model that generates 
$ is not unique. Since $ is the same for both impu- 
rity and lattice model, they have the same diagrammatic 
expansion. We therefore embed a single unit cell (a "o 
/ molecule") as impurity (H\ oc ) in an effective medium 
(i? mc( j) which will be determined self-consistently: 



-f^imp — H-loc "b Hraed 



(12) 



cr a 
a 

i? m ed = E^ a L« feCT + E(^4« feCT +H.c.) . (14) 



ft (7 



fc<T 



(Formally one can consider the action of the lattice 
model and integrate out the non-local part to arrive at 
an impurity action which is afterwards modelled by an 
Hamiltonian.Ej) Note that this choice for the impurity 
model differs qualitatively from tho-u suaL one for the pe- 
riodic Anderson model at U c = OEjESEZI In the case of 
free conduction electrons, only a single self energy ex- 
ists (for the /-electrons) and, hence, only the /-orbital 
is coupled to an effective medium. When the conduc- 
tion electrons are correlated they have a self energy, as 
well. Therefore, we include a c-orbital in the local part 
of -ffimp. We need only a single effective medium, i? mo d, 
although there are two Hubbard interactions in the origi- 
nal Hamiltonian (Q) . This is due to the absence of direct 
/-/ hopping: The electrons explore their environment 
only via the c-orbitals. The two interactions merely show 
up in the internal structure of the impurity which con- 
sists of two orbitals, one of which couples to the medium. 
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This situation is similar to the one encountered in the ex- 
tended Hubbard model in Ref. 

We show now that the self-consistency equation ( |ll| ) 
can be fulfilled by our impurity model (|12j) . The impurity 
Green's functions 



are given by 



Gf(z) Qf c {z) 
Gcf{z) g c (z) 



-V-^cf(z) 



u-ef-Ef(z) 
-V-T,f c (z) uj - e c - A(z) - T, c (z) 



(15) 



Q{z) = 1 , 



where 



*(*) = E 



\W k \ 2 

z - Eh 



(16) 



(17) 



Equating Q to the local Green's function G\ oc of the lat- 
tice model (f[o]), which one obtains analytically due to 
the density of states (||), and using that both, A(z) and 
G(z) ~ l/z at large \z\, we find: 



i c = e c - p 
A(z) = -Gi oc ,cO) 



(18) 
(19) 

(20) 



The medium is determined by the c-Green's function 
only, reflecting that there is no direct /-/ hopping. 

The parameters in the impurity model ( pL2|) is thus 
fixed, in particular the medium A(z) is determined by 
the solution of the impurity model Q c {z) = Gi oc c (z). 



C. Solving the impurity model 

We solve the impurity model (jlfl^pby extending the 
non-crossing approximation (NGAjOLj to the case of 
more than two ionic propagators J 3 ir 3 ! This approach has 
been applied successfully to the finite-?/ impurity Ander- 
son mode£3 where it has been shown that neglecting ver- 
tex corrections slightly underestimates the Kondo tem- 
perature, and to lie Emery model within the dynamical 
mean- field theoryEj 

Denoting the eigenstates of the local part, H\ oc by 
|m) (with n m — 0...4 particles), the impurity Hamil- 
tonian ([[2]) is expressed in terms of Hubbard operators 
X mn = \m) (n\ as 
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m—l 

-Hmed = E Eka\ a QLk<j 

E (WkU^ n!7 X mn a ka + H.c. 



(21) 



(22) 



with U^ na = (m|4|n)._ 

For each state \m) a ionic propagator is introduced 



R m {z) 



1 



z E m S m {z) 



(23) 



with spectral density p m (x) — — Im R m (z + i0 + )/7r. We 
assume that the corresponding self-energies S are diago- 
nal in the local basis and evaluate them in self-consistent 
perturbation theory to second order in the hybridization 
W as in the usual NCA: 



5 m (z)=E(|C/ I c „„J 2 + |C/^ m J 2 )x 



de f{r/ mn e)A(e)R n (z + r] mn e) . (24) 



Here rj mn = — 1(+1) if the particle number in |m) is 
higher (lower) than in \n), f(z) — [exp(/3z) + 1] _1 is 
the Fermi function, and A = — Im A/tt. The c- and f- 
Green's functions are given by 



G c (f){z) — \ U^^\ 2 ((X mn ; X nm )) 2 



(25) 



Within the NCA the Green's furictioHSi {{X mn ; X nm )) 2 
are expressed by the ionic propagators" 



( {X mn , X nm } ) z — 
1 



(26) 



dx e 0x {p m {x)R n {x + z) -p n (x)Rm(x- z)} (27) 



and 



Z = Y,J dx e-^p m (x) 



(28) 



In the symmetric, half-filled case there arc only 6 inde- 
pendent propagators due to particlc-holc and spin sym- 
metry. The coupled integral equations ( |23] ) and (|^) axe 
solved numerically by introducing defect propagators! 3 ^ 
and making use of the fast Fourier transformation^ 



III. RESULTS 

In the following we consider the symmetric model (e c = 
-C4/2, e/ = -Uf/2) at half-filling (n c + n f = 2). Due 
to particle-hole symmetry the chemical potential is 0. 

We chose Uf = 5 in all our calculations so that the 
/-level is well outside of the conduction band. Our in- 
vestigations were restricted mostly to those values of U c 
for which the Hubbard model for the conduction elec- 
trons is metallic (for the semielliptic density of states (||) 
we found U CT it ~ 1.8). The reason is that in deriving 
the self-consistency equations for the impurity model we 
assumed a paramagnetic state. Thus we typically chose 
U c = 0.5 . . . 2.0 and V = 0.1... 0.4. These values lead to 
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V 




J K 


0.1 


1.3 x 10" 44 


2.7 x 10"^ 


0.2 


2.5 x 10~ 12 


5.3 x 10~ 5 


0.3 


3.1 x 10~ 6 


7.3 x 10~ 4 


0.4 


4.9 x 10 -4 


1.1 x 10~ 2 



TABLE I. The estimated Kondo temperatures for uncor- 
related conduction electnme. T™ p refers to the symmetric 
impurity Anderson modelJSZI Tjf* includes the lattice enliance- 
ment factor of 2 in the exponent for the lattice case£3'C3 



small exchange couplings J ~ 0.01. ..0.2 between two- 
particle singlet and triplet state of the c-f molecule (see 
below) . We deliberately chose these small values for V in 
order to obtain "Kondo temperatures" which are small 
compared to the bare band-width, see Table |. 



A. Single-particle spectrum of the molecule 



They correspond to the lower Hubbard bands of the c- 
and /-subsystem which are shifted by the hybridization 
V. 

The first excited state is the two-electron triplet state 
(spin excitation) 



1 



l*T> = ^(^+cI/f)|0) 
with excitation energy 



(34) 



AE = E T -E s = ^ U f + U ^ + W*- U f + U c 
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8V 2 



u f + u c 



(35) 



At the temperatures that we investigate in the following 
(T = 0.5 • • • 0.001) only these two states contribute as 
initial states, \m). The resulting photoemission spectrum 
of the molecule thus consists of two double peaks at z c , 
z c + AE, and Zf, Zf + AE and the weight of the peaks 
shifted by AE goes to as T -> 0. 



The following discussions focus on the one-particle 
spectral function. It is instructive to investigate this 
quantity first for the local problem given by H\ oc in 
Eq. ( p"2| ) in the symmetric case. Due to symmetry, we 
consider the photoemission spectrum only. Denoting the 
eigenstates and eigenenergies of -ffi oc by |m) and E m , the 
one-particle Green's function and spectral function are 
given by (fi = 0) 



G c (z) 



Z 



e -? E "> x 



mn 

m\cf\n) (n\c^\m) 



(m\c\\ri) (n|c-|-|m) 



A c (z) 



z-{E n -E m ) z+(E n -E m ) 
-Im G c (z + i0+) . 



(29) 



(30) 



At zero temperature only the ground state contributes 
to the sum over m. It contains two electrons forming a 
singlet 



|4>c 



Ea = 



— (c f f f - c f f f ) + — IcV - f f f f ) 



u f + u c 4 l(U f + U c ) 
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|0) (31) 



(32) 



The photoemission spectrum is obtained by removing a 
particle, hence the final states are bonding and antibond- 
ing combination of the c- and /-orbital. To lowest order 
in V the transition energies are given by 



Uc 

2 



3U f - 5U C 



(U c + Uf)(U f -U c ) 



V 2 



_ U f 5U f -3U c 
f 2 [U c +Uf)(U f -U B ) 



v z 



(33) 



B. Single-particle spectrum of the 
Anderson-Hubbard model 

We first consider the spectral density of the conduction 
electrons 

1 



A c (z) = Im G loc , c (z + i0+) 



-Vim G c (k,z + i0+) 

IT ' 



(36) 



which corresponds to photoemission and inverse photoe- 
mission spectra. A typical result for A c (z) is shown in 
Fig. |l] for different temperatures (U c = 1, V = 0.2). At 
high temperatures (T = 0.5) we obtain two broad max- 
ima located at ~ ±[/ c /2 which correspond to upper and 
lower Hubbard band of the c-subsystem. When lower- 
ing the temperature (T < 0.1), a peak at the chemical 
potential (z — 0) arises. This is the well-known quasipar- 
ticle peak belpagkia to the Hubbard model of the host 
electrons, H c P^ PtI 25 ! In this temperature regime the /- 
and c-subsystems are almost independent. The influence 
of the /-states on the c-spectra is given only by the tiny 
structure at z ~ ±Uf/2 — ±2.5. The separation of the 
two subsystems is indicated in Fig. ^ as well where we 
compare the spectral weight at the chemical potential for 
different values of the c-f hybridization V and the pure 
c- Hubbard model (V = 0): At high to moderate tem- 
peratures the spectral weight of the Anderson-Hubbard 
model follows the one of the pure Hubbard model. 

In Fig. PJ we also see that this behavior does not ex- 
tend to low temperatures. At a certain temperature, 
which depends on V, the spectral weight no longer fol- 
lows the quasiparticle peak of the Hubbard model but 
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FIG. 1. Spectral density A c {z) for U C = 1,V = 0.2, U f = 5 
at different temperatures 




FIG. 2. Spectral weight of the conduction electrons at the 
chemical potential, A c (0), for U c = 1 vs. temperature 



drops to zero. As seen in Fig. [j] indeed a gap occurs and 
sharp structures emerge close to the chemical potential 
when the temperature falls below T < 0.04. This resem- 
bles the Anderson model with uncorrelated conduction 
electrons.E3 There, at T < Tq where To is a characteris- 
tic temperature related to the Kondo temperature, the 
Kondo effect leads to a resonance at the chemical po- 
tential. These dynamically generated local states cross 
the conduction-band states and one finds a splitting of 
the conduction band with a gap at the position of the 
resonance. Due to particle-hole symmetry, this feature 
occurs at the-chemical potential and the system becomes 
an insulator SB 

A similar behavior is found in the case of interacting 
conduction electrons. In order to see that indeed flat 
bands occur close to the chemical potential we inspect 
the fc-dependent spectral function 



A c (k,z) = Im G e (k,z) 

7T 



(37) 



From the self-consistency equation (|Tl]) and Eq. (Q) we 
find 



G c (k,z) 



1 



G iocJ z ) ~ e k + iG loc , c (z) 



(38) 



Since A c (k, z) depends on k only via we plot A c (k, z) 
for different e& = —1 ... 1 in Fig. |[ At high temperatures 
(T = 0.3, Fig. |a.) the features are very broad and we 
basically see the Hubbard bands for the c-subsystem at 
~ U c /2 as well as for the /-subsystem at ~ Uf/2 which 
weakly admix to the c spectra. At intermediate temper- 
atures (T = 0.1 . . . 0.02, T = 0.05 is shown in Fig. §6.) 
spectral weight is found at the chemical potential as well 
and we may trace the quasiparticle band of the Hubbard 
model. When the temperature is decreased below 0.015 
(T = 0.005, Fig. ||c.) these peaks split, a small gap opens 
and the newly emerged peaks show a weak dispersion 
at the chemical potential. The transition occurs rather 
quickly: Whereas at T = 0.015 only the peak at e& = is 
split and the band follows the quasiparticle band of the 
Hubbard band elsewhere, two separated bands already 
emerged at T = 0.01. Due to their weak dispersion at 
the chemical potential we expect that they will lead to 
heavy bands upon doping. At higher energies these new 
bands merge in the previous quasiparticle bands of the 
Hubbard model. 

These findings fit qualitatively to the scenario of the 
Anderson modeL-svith uncorrelated conduction electrons 
described aboveEa and-jto the results of the slaved-boson 
mean-field treatmentJi!J In the latter, the free conduction 
bands hybridize with a strongly renormalized (reduced) 
coupling to an effective /-level at the chemical poten- 
tial. This leads to weakly dispersive bands at low ener- 
gies merging in the original bands at high energies. In 
the Anderson-Hubbard model the new bands merge in 
the former quasiparticle band of the c-Hubbard model at 
high energies. We conclude that these quasiparticles take 
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FIG. 3. Spectral function A c (k, z) for U c = 1, V = 0.2 and 
different values of e(k). a. T = 0.3, b. T = 0.05, c. T = 0.005 
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FIG. 4. Spectral density A c {z) for U c = 1.5, V = 0.2, 
Uf = 5 at different temperatures 

the role of the free conduction electrons in the dynamical 
screening of the /-moments leading to the resonance at 
the chemical potential. This interpretation does not con- 
tradict previous results on a single impurity embedded 
in a correlated host.EJ There it turned out that,a vari- 
ational ansatz in the spirit of Varma and Yafetcil which 
uses quasiparticles for the screening of the /-moment in- 
stead of bare electrons is not sufficient to find the correct 
Kondo temperature. This holds due to importance of the 
renormalization of the c-f exchange interaction and does 
not imply that the quasiparticle picture is not valid in 
describing the screening process. 

When the interaction strength of the conduction elec- 
trons is increased to U c — 1.5, no qualitative changes 
occur at first sight. In Fig. [| we show the c-spectra at 
various temperatures for U c — 1.5 and V = 0.2. Again, 
the c- and /-subsystem are separated at high tempera- 
tures and a quasiparticle peak of the c-Hubbard model 
evolves first. At low temperatures a Kondo resonance 
is formed at the chemical potential with hybridizes with 
the quasiparticle band, a gap opens and we find bands 
with a weak dispersion. When compared to U c = 1, we 
find that the gap has increased. This indicates that the 
hybridization of the quasiparticle band with the dynam- 
ically generated states has increased. 

However, as demonstrated in Fig. |^, a quasiparticle 
peak does not always occur in the intermediate temper- 
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FIG. 5. Spectral density A c (z) for U c = 1.5, V = 0.3, 
Uf = 5 at different temperatures 



ature range even though the bare conduction electrons 
are metallic: Increasing the c-f mixing to V — 0.3, the 
system becomes directly insulating although U c — 1.5 < 
U cr if Figure || illustrates that it depends on the value 
of V whether the quasiparticle peak shows up or not. 
We conclude that part of the strong correlation on the 
/-orbital is effectively inherited by the c-orbital via the 
hybridization V. A similar effect has been observed and 
discussed for a different model in Ref. |42| 

Although no quasiparticle peak corresponding to the 
conduction electrons emerges in Fig. |^, we still recover 
the Anderson scenario described above when decreasing 
the temperature, i.e., a gap opens and peaks arise close 
to it at low temperatures T ~ 0.04. It appears that we 
do not need pre-formed quasiparticles in order to screen 
the /-moments which in turn leads to the Kondo reso- 
nance. Note however, that there is finite spectral weight 
at the chemical potential. This is better seen in the k- 
resolved spectral function A c (k, z) in Fig. ^ where we 
display A c (k,z) for two different temperatures for the 
Anderson-Hubbard (U c = 1.5, V = 0.3) and the Hubbard 
model (U = 1.5). For the Hubbard model we find a peak 
crossing the chemical potential at T = 0.07 (Fig. 0a.) . 
It corresponds to the quasiparticle peak of the Hubbard 
model being formed at this temperature. The Anderson- 
Hubbard model, however, exhibits no structure crossing 
the chemical potential. Note that upper and lower Hub- 



FIG. 6. Spectral weight of the conduction electrons at the 
chemical potential, A c (0), for U c = 1.5 vs. temperature 



bard bands of the Anderson- Hubbard model are in agree- 
ment with those of the pure Hubbard model. At T = 0.04 
(Fig. 0o.) the Anderson-Hubbard model shows two flat 
bands at the chemical potential. The peaks are most 
pronounced close to chemical potential. In contrast to 
the case of small U c (U c — 1 discussed above) the corre- 
sponding bands do not merge the quasiparticle band of 
the Hubbard model. This indicates a large effective hy- 
bridization between quasiparticle and dynamically gener- 
ated states at the chemical potential. To a lesser degree 
this behavior is also observed for U c = 1.5 and V = 0.2 
(A c (k, z) not shown). We conclude that the effective hy- 
bridization increases as both U c and V increase. Consid- 
ering again Fig. 0a., we find that the difference between 
Hubbard and Anderson-Hubbard model at high tempera- 
tures is related to the onset of the heavy bands which are 
thus already seen at comparatively high temperatures. 

We finally turn to the case where the bare c-subsystem 
is insulating, i.e., U c > C/ C rit- For U c = 2 a shoulder 
in the spectral function emerges at the edge Hubbard 
band towards the chemical potential when the temper- 
ature is lowered (Fig. ||). That this feature is indeed 
caused by the /-subsystem is demonstrated in Fig. ^ 
where we compare A c {k,z) for the Anderson-Hubbard 
model and the pure Hubbard model. Whereas the struc- 
tures corresponding to the Hubbard bands roughly agree, 
the Anderson-Hubbard model shows an additional peak 
at the low-energy edge of the Hubbard bands. This is un- 
expected because the c-subsystem is insulating and pro- 
vides no quasiparticles which could screen the /-moment 
and the resulting spectra shoul d rese mble the one of the 
c-f molecule described in Sec. Ill A . We do not believe 
that this shoulder is spurious since increasing the energy 
resolution which allows to proceed to lower temperatures 
did not change the shoulder. However, we can not ex- 
clude a principle failure of the extended NCA as it is 
well-known that the NCA-fails to converge when the sys- 
tem becomes insulating.CJ On the other hand, the spec- 
tral weight at the chemical potential is not exactly zero. 
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FIG. 7. Comparison of the spectral function A c (k, z) of the 
Anderson-Hubbard model (U c = 1.5, V = 0.2, solid lines) to 
the Hubbard model (U c = 1.5, dashed lines, multiplied by 
-1) at a. T = 0.07 and b. T = 0.04. 



FIG. 8. Spectral density A c (z) for U c = 2, V = 0.3, U f = 5 



Similar to the case U c = 1.5, V = 0.3 discussed above, 
this small, but finite spectral weight could be sufficient 
to dynamically generate states at the chemical poten- 
tial which result in a heavy band via hybridization as 
previously. To clarify this situation, it is mandatory to 
employ other numerical methods for solving the impurity 
model. If the picture presented above is valid, the effec- 
tive hybridization must be large compared to the metallic 
cases so that the new bands are pushed towards the Hub- 
bard bands. Note that the new peaks approach, but not 
merge in the Hubbard bands whereas they merged in the 
quasiparticle bands in the metallic cases (U c — 1, and 
U c = 1.5, V < 0.2). 



C. Hybridization gap 

In this section we focus on the gap which opens at low 
temperatures. We determine the size of the gap, A, as 
twice the distance between zero frequency.-and the first 
maximum close to the chemical potential. c3 In Fig. [l(], 
where we plot the gap A vs. temperature for different 
hybridizations V, we observe that the gap opens quickly 
when the temperature is reduced and we define Tq as the 
temperature where the gap opens. One expects that Tq 
is related to the energy difference, AE, of singlet and 
triplet states in the local problem, see Sec. Ill A. In- 
deed, we observe in Figs. || and || that the gap opens 




FIG. 9. Comparison of the c spectral function A c (k,z) of 
the Anderson-Hubbard model (U c = 2, V — 0.3, Uf = 5, solid 
lines) to the Hubbard model (U c = 2, dashed lines, multiplied 
by -1) at T = 0.065 




0.20 



FIG. 10. Gap in the spectral density vs. temperature. The 
gaps are obtained from the two maxima in A c (z) for U c = 1 
and different hybridizations V 
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FIG. 11. Splitting of the / peak in the molecule and lattice, 
and T vs. U c at V = 0.3 



roughly in the same temperature regime where the /- 
peak at z ~ ±Z7//2 splits, i.e., where the states |$g) and 
|\&t) become distinguishable in the photoemission of the 
lattice. We extract this splitting of the /-peak, AE\ at 
from the spectral function of the lattice model (if visible) 
and compare it to the corresponding splitting, A_Ei oc for 
the c — / molecule [cf. Eq. (pq)] in Fig. [TT] for different 
values of U c and V = 0.3. Surprisingly we find that while 
both splittings are of the same order of magnitude, they 
depend differently on U c : When U c increases, A£i at in- 
creases, whereas AEi oc decreases. 

In Fig. [12] we plot the gap vs. temperature for different 
values of U c and V = 0.2, 0.3. In general, the tempera- 
ture To at which the gap sets in increases with U c , but 
deviations occur: When V — 0.3, U c = 1.5 and 1.7 do 
not fit into this scheme. However, one should bear in 
mind that the onset of the gap formation as we measure 
it, depends also on the shape of the spectrum. As is seen 
from Fig. [H], To is of the same order of magnitude as 
Ai?i a t and its general behavior agrees with AE\ a t and is 
thus opposed to AE\ oc . 

We now turn to the size of the gap. It has been shown 
for the Anderson model with free conduction electrons 
that this quantity determines the low-temperatur* 
modynamics in the limit of infinite dimensions.! 
The gaps in the local spin and charge excitation spec- 
tra were found to be of the same order of magnitude as 
the gap in the density of states. The latter thus provides 
a measure of the low-temperature scale ( "Kondo temper- 
ature"). From Fig. |lj we extract that the size of the gap 
has not yet converged at the lowest temperatures where 
our extended NCA ceases to converge. Nevertheless, we 
can draw some qualitative conclusions: The size of the 
gap increases systematically as U c increases. This is also 
deduced from Fig. 13 where we plotted the gap at the 



lowest temperature we could reach for each pair (U c , V) 
vs. the strength U c of the Hubbard interaction of the con- 
duction electrons. This procedure implies that the points 
shown correspond to different temperatures. Note that 







the magnitude of A varies much stronger with U c and V 
compared to Tq or AE. 



< 




0.04 



0.00 0.02 0.04 0.06 0.08 

T 

FIG. 12. Gap in the spectral density vs. temperature. The 
gaps are obtained from the two maxima in A c (z) for different 
interaction strength U c and a. V = 0.2 and b. V — 0.3 




FIG. 13. Gap in the spectral density vs. U c - The gaps 
correspond to the lowest temperatures reached. 



IV. CONCLUSIONS 

In conclusion, we studied the influence of interactions 
among the conduction electrons on the low-temperature 
behavior of the periodic Anderson model. In the dynam- 
ical mean-field theory the model is mapped onto a gener- 
alized Anderson impurity model that couples the orbitals 
of a single unit cell to an effective medium which has to 
be determined self consistently. The impurity model was 
solved numerically by an extended non-crossing approx- 
imation. 

For weakly interacting conduction electrons we found 
that at high temperatures the c- and /-subsystems are al- 
most separated as in the case of free conduction electrons. 
Decreasing temperature then first leads to the formation 
of quasiparticles in the c-subsystem as in the bare Hub- 
bard model. When the temperature is reduced further, 
the quasiparticle band splits, a tiny gap opens and the 
system turns into an insulator. As in the case of free con- 
duction electrons, the gap is formed by the level crossing 
of the quasiparticle (Hubbard) band and the resonance 
which arises at the chemical potential from the Kondo- 
like screening of the /-moments. We observed that the 
quasiparticles play an essential role in this screening. The 
resulting two bands have a weak dispersion close to the 
gap and we expect that they turn into heavy bands upon 
doping. 

When the correlations of the conduction electrons be- 
come stronger, the low-temperature gap increases. This 
can be interpreted as increasing the effective hybridiza- 
tion between the quasiparticle states at the chemical po- 
tential and the dynamically generated states which leads 
to a larger gap and is qualitatively in agreement with re- 
sults found for impurity models. For the latter case, the 
main effect of the (small) interaction was to renormalize 
and increase the exchange interaction. til 

It turned out, however, that pre-formed quasiparticles 
within the c-subsystem are not prerequisite for the emer- 
gence of heavy bands. When increasing the c-f hybridiza- 
tion, the c orbitals seem to inherit correlations from the 
/ orbitals and a quasiparticle peak is no longer formed 
in the spectral density at intermediate temperatures, al- 
though the spectral weight at the chemical potential does 
not vanish. Nevertheless, we observed heavy bands at low 
temperatures in the one-particle spectra. 

Even when the (bare) c-system is insulating, it is influ- 
enced by the /-system at low temperatures: A shoulder 
forms at the edge of the Hubbard band which shows weak 
dispersion. This is surprising since the bare c-system pro- 
vides no quasiparticles which could screen the /-moments 
and the c spectral weight close to the chemical potential 
is small. However, we can not exclude that this result is 
an artefact of our method to solve the impurity problem 
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and further investigations are necessary to decide upon 
this question. 

Finally, we investigated the hybridization gap which 
occurs at low temperatures. When increasing the 
Coulomb interaction among the conduction electrons the 
temperature at which the gap opens increases. This tem- 
perature is related to the splitting of the /-peak in the 
spectrum which results from singlet and triplet states 
in the c-f molecule. However, we found that this split- 
ting depends oppositely on the interaction strength in the 
lattice and in the molecule. The low-temperature ther- 
modynamics in infinite dimensions scales with the size of 
the gap.Ea This quantity therefore provides a measure for 
the "Kondo temperature." In agreement with impurity 
models, this gap increases as the correlations among the 
conduction electrons become stronger. 
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